第 38 卷 第 8 期 工程 热 物理 学 报 ‘Vol.38, No.8 
2017 年 8 月 JOURNAL OF ENGINEERING THERMOPHYSICS Aug., 2017 


心静 平 设计 与 变 工 况 性 能 研究 


宋 艳 蔷 "” 谭 金 ” 黄 典 贵 


(1. 上 海 理工 大 学 能 源 与 动力 工程 学 院 ， 上 海 200093; 
2. 上 海 市 动力 工程 多 相 流动 与 传 热 重 点 实验 室 ， 上 海 “200093) 


摘 要 本文 以 R123 为 流动 工 质 ， 进 行 离心 透 平 级 的 一 维 热力 计算 ， 并 通过 数值 模拟 分 析 ORC 离心 透 平 级 内 流动 情 
况 及 变 工 况 性 能 . 结果 表明 ;在 设计 工 况 下 ,数值 模 拟 结果 与 一 维 设计 结果 基本 一 致 功率 与 效率 偏差 在 1 % 以 内 ， 符 合 
设计 要 求 . 在 变 工 况 下 ， 进 汽 压力 在 0.33~0.88 MPa 变化 时 ,效率 随 进 汽 压力 增 大 先 增 大 后 减 小 ,流量 和 功率 随 进 汽 压 
力 增 大 而 增 大 ; 背 压 在 0.10~0.18 MPa 变化 , 在 临界 工 况 下 运行 , 效率 和 功率 随 着 青 压 的 增 大 而 减 小 , 流量 保持 不 变 ; 初 
温 改变 对 透 平 性 能 影响 较 小 , 效率 变化 范围 在 1 % 以 内 ; 转速 变化 时 , 流量 不 变 , 效率 先 增 大 后 减 小 , 设计 转速 下 效率 最 高. 
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Abstract 
by using R123 as working fluid. Combining with CFD technology, the aerodynamic performance 


In this paper, the one-dimensional design of the centrifugal turbine was conducted 


analysis of the centrifugal turbine stage was carried out both in design and off-design conditions. 
The results show: in the design condition, the CFD results mostly agree with the one-dimensional 
design; relative errors of power and efficiency are within 1%, efficiency and power of the turbine 
satisfy the design requirements. In the off-design conditions, when the inlet pressure increased 
from 0.33 MPa to 0.88 MPa, the efficiency increased at first, then decreased; mass flow and power 
increased with the increasing of the inlet pressure. When the back pressure increased from 0.10 MPa 
to 0.18 MPa which is smaller than the critical back pressure, efficiency and power decrease with the 
The variation of inlet temperature 
When the rotor 


speed varies, efficiency increases with the increasing of the rotor speed, reaches a maximum value at 


increasing of the back pressure, mass flow maintains constant. 


showed less effect on stage performance, the variation of efficiency is within 1%. 


design speed, and then decreases with any further increase in rotor speed, mass flow maintains constant. 
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在 我 国 , 工业 余热 、 太 阳 热 能 、 地 热能 等 低 品位 
热源 强 藏 量 大 ， 分布 广 ， 低 品位 热源 科学 利用 对 能 
源 开发 具有 重大 的 意义 . 但 由 于 常用 的 工 质 水 薰 气 
和 理想 气体 等 沸点 较 高 ， 低 温 热源 进行 热切 转换 效 
率 低 . 近年 来 , 由 于 有 机 工 质 低 沸 点 的 特性 , 有 机 朗 
肯 循 环 (ORC) 成 为 低温 热源 利用 研究 的 热点 。 

透 平 是 ORC 热 功 转 换 的 关键 设备 , 是 系统 
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高 效 运 行 的 保证 。 有 机 工 质 分 子 量 大 ， 沸 点 低 ， 体 
积 变化 大 , 声速 小 等 不 同 于 理想 气体 的 特殊 热 物性 ， 
使 透 平 内 流动 也 更 为 复杂 ， ORC 透 平 的 研 
究 显 得 愈加 重要 . 刘 广 彬 日 采用 R123 作为 循环 工 
质 , 以 涡 旋 膨 胀 机 作为 能 量 回 收 机 械 , 建立 了 小 型 低 
温 余 热 发 电 试验 系统 , 分 析 了 相同 热源 不 同 流量 的 
两 种 工 况 下 ， 膨胀 机 转速 对 系统 性 能 的 影响 .Wang 
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W, Wu Y T Ps3l 研究 单 螺杆 膨胀 机 用 于 余热 回收 
发 电 , 最 低 可 利用 余热 温度 可 低 至 80"C. 薄 泽 民 欠 
设计 了 150 kW 向 心 透 平 ， 采 用 NUMECA 软件 进 
行 数值 模拟 ,并 分 析 了 其 变 工 况 性 能 。 李 艳 [5.9 设 
计 了 高 膨胀 比 有 机 工 质 向 心 透 平 ， 进 行 了 数值 模拟 
和 实验 .意大利 米兰 理工 大 学 的 G.Persico，M.pini 
等 中 研究 了 离心 透 平 流动 规律 及 应 用 范围 ,对 离心 
透 平 叶 栅 的 气动 性 能 进行 了 数值 模拟 研究 

本 课题 组 李 银 各 四 以 完全 气体 为 工 质 , 研究 了 
一 维 级 效率 的 影响 因素 ， 为 离心 透 平 的 设计 提供 了 
依据 . 谭 从 外 对 以 完全 气体 为 工 质 的 离心 透 平 进行 
了 气动 设计 和 优化 ， 并 采用 数值 模拟 方法 ， 对 所 设 
计 的 离心 透 平 级 进行 设计 工 况 及 变 工 况 三 维稳 态 流 
场 分 析 . 黄 典 贵 、 谭 从 [9 设计 了 一 种 新 型 的 离心 透 
平 , 结构 如 图 1 所 示 . 气流 从 内 部 环形 径 向 进 汽 , 脱 
胀 做 功 后 沿 外 径 径 向 流出 ， 通 流 截面 沿 着 流向 膨胀 
过 程 逐渐 增 大 。 


图 1 离心 式 透 平 结 构 示 意图 


Fig. 1 The schematic of centrifugal turbine 


有 机 朋 肯 循环 采用 该 离心 透 平 与 传统 的 向 心 透 
平 相 比 有 两 点 优势 : 1、 有 机 工 质 在 膨胀 时 ， 相 对 于 
理想 气体 比 容 变化 大 (该 设计 几何 与 气动 相 匹配 , 叶 
高 变化 小 ) 2、 有 机 工 质 分 子 量 大 ,声速 低 ， 容 易 形 
成 超 音 速 流动 ， 引 起 激 波 和 气流 偏转 。 而 离心 透 平 
可 采用 多 级 设计 ， 各 级 烩 降 合理 分 配 ， 避免 在 马 替 
数 过 高 的 工 况 下 运行 . 同时 多 级 可 利用 余 速 和 重 热 ， 
提高 整 机 热效率 . 本 文 以 R123 作为 流动 工 质 , 对 大 
膨胀 比 离心 透 平 进行 气动 设计 ， 并 数值 模拟 研究 其 
流动 特性 及 变 工 况 性 能 。 


1 离心 透 平 气动 设计 
根据 ORC 离心 透 平 的 工 质 特性 和 透 平 流动 特 
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征 , 借鉴 轴 流 和 向 心 透 平 气动 设计 方法 044234， 建立 
流动 损失 模型 , 以 轮 周 效率 最 大 为 目标 , 热 物性 参数 
从 实际 气体 物性 软件 REFPROP9.0 调用 , 编写 一 维 
热力 设计 程序 .本文 以 某 太阳 能 热源 有 机 朗 肯 循环 
向 心 透 平 热力 参数 见 表 1， 以 R123 为 流动 工 质 , 根 
据 编制 的 一 维 热力 计算 程序 , 设计 离心 透 平 , 以 轮 周 
效率 最 大 为 目标 ,采用 筛选 法 ,确定 设计 几何 参数 ， 
如 表 2 所 示 ， 


表 1 离心 透 平 热力 参数 
Table 1 The parameters of the centrifugal 


turbine 
热力 参数 设计 值 
进口 总 温 To0/K 373.23 
进口 总 压 Po/MPa 0.78564 
出 口 背 压 P/MPa 0.11 
设计 转速 n/r-min-! 10000 
质量 流量 G/kgs-! 12.4 
透 平 功率 W /kW 300 
轮 周 效率 / 吧 84 


表 2 离心 透 平 的 主要 几何 参数 
Table 2 The centrifugal turbine geometry 


热力 参数 设计 值 
静 叶 入 口 直径 Di /mm 278 
静 叶 出 口 直径 D2/mm 310 
动 叶 入 口 直径 D3/mm 314 
动 叶 出 口 直 径 D4 /mm 346 
叶 高 H/mm 21.7 
进口 气流 角 a/(?) 12 


2 离心 透 平 的 数值 模拟 

根据 离心 透 平 的 一 维 热力 设计 参数 ， 采用 
Ansys-Workbench 工作 平台 系列 软件 ， 分 别 对 动静 
叶 进 行 叶 型 初步 造型 设计 、 叶 型 参数 化 、 叶 栅 流 道 
网 格 划 分 、 数 值 模拟 计算 ， 在 初步 模拟 的 结果 基础 
上 采用 梯度 算法 对 叶 型 进行 优化 ,循环 数值 模拟 计 
算 之 至 功率 与 效率 达到 设计 要 求 ， 然 后 进行 变 工 况 
性 能 分 析 . 具体 流程 如 图 2 所 示 。 
2.1 叶片 造型 

离心 透 平 通 流 截面 和 工 质 比 容 沿 流向 方向 同时 
增 大 ， 因 此 叶片 可 采用 等 叶 高 直 叶 片 ， 但 由 于 环 列 
叶 栅 流 道 是 渐 扩 的 ， 叶 型 无 法 采用 传统 的 轴 流 叶 型 ， 
设计 满足 离心 透 平 特性 的 气动 经 型 是 一 个 重要 的 工 
作 . 本文 在 Ansys 软件 中 BladeGen 模块 中 设计 叶 
型 .以 静 叶 为 例 予 以 简要 说 明 : 叶片 为 等 叶 高 直 叶 
片 ， 根 据 表 2 确定 叶片 的 进口 径 向 位 置 与 叶 高 确 
定 叶 型 子午 面 。 叶 片 截面 由 前 缘 、 尾 缘 、 叶 贫 和 叶 
背 四 条 曲线 组 成 ， 前 缘 和 后 缘 采 用 圆 弧 ， 叶 背 、 叶 
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倪 则 由 中 弧 线 和合 加 厚度 的 方法 来 构造 。 如 图 3 所 
示 ， 中 弧 线 是 采用 3 次 Bezier 曲线 来 表示 ， 四 个 
点 代表 不 同 径 向 位 置 的 中 弧 线 与 切 向 的 夹 角 ， 横 坐 
标 代表 的 是 径 向 相对 位 置 ， 纵 坐标 表示 相应 的 几何 
角 , 通过 改变 中 间 2 个 控制 点 的 坐标 可 改变 中 弧 线 
形状 。 沿 中 弧 线 的 厚度 分 布 也 采用 3 次 Bezier 曲 
线 来 表达 ， 如 图 4 所 示 ， 前 缘 和 尾 缘 厚 度 为 圆 弧 直 
径 ， 经 过 经 验 值 确定 ， 通 过 改变 中 间 2 个 控制 点 
的 坐标 来 得 到 中 弧 线 不 同 径 向 相对 位 置 的 厚度 值 . 
最 终 由 样 条 曲线 光滑 连接 生成 叶 贫 、 叶 背 曲 线 . 同 
理 生 成 动 叶 叶 型 ， 图 5 所 示 为 动静 叶 三 维 几 何 结 
构图 。 
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各 用 梯度 算法 优化 


达到 设计 要 求 


图 2 离心 透 平 的 叶 型 设计 和 优化 流程 图 


Fig. 2 Numerical simulation of centrifugal turbine 
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图 3 中 弧 线 切线 角 分 布 
Fig. 3 Mean line tangential angular distribution of 
radial chord 


2.2 网 格 划分 与 边界 条 件 的 设置 

根据 叶 栅 的 流动 特性 ， 在 数值 模拟 时 选取 单个 
流 道 为 研究 对 象 简化 计算 ， 周期 性 边界 条 件 ， 近 壁 
面 为 无 滑 移 固 体 壁 面 , 网 格 采 用 ANSYS-TurboGrid 
生成 结构 网 格 . 满 流 模型 选择 k-e 双方 程 模型 , 标准 


壁面 函数 . 工 质 物性 选用 实际 气体 立方 型 S-R-K 状 
态 方 程 计算 ,进口 边界 条 件 给 定 总 温和 总 压 ， 出 口 
边界 给 定 静 压 , 设计 转速 为 10000 r/min。 
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图 4 叶片 厚度 分 布 
Fig. 4 _ Thickness distribution of radial chord 


图 5 叶片 几何 模型 


Fig. 5 3-D model of stator and rotor 


2.3 数值 模拟 结果 

根据 数值 模拟 结果 是 否 达到 设计 效率 及 功率 ， 
采用 梯度 优化 算法 ,通过 调整 叶片 的 中 弧 线 切线 角 
和 叶片 厚度 沿 径 向 分 布 ， 优 化 叶 型 ， 数 值 模拟 重复 
计算 直至 达到 设计 要 求 。 

表 3 给 出 了 所 设计 的 离心 透 平一 维 设计 与 CFD 
数值 模拟 的 结果 , 在 设计 工 况 下 , 热力 计算 和 数值 模 
拟 结果 误差 在 1% 以 内 ， 切 率 和 效率 达到 设计 要 求 ， 
符合 预期 。 


表 3 设计 值 与 模拟 结果 误差 表 


Table 3 The contrast of design value and 


CFD value 
变量 设计 值 模拟 值 
G/kg:s-! 12.4 12.55 
W/kW 344 343 
人 2 0.6 0.62 
Xa 0.63 0.6 


n/% 84 83 
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图 5 所 示 , 级 的 动静 叶片 均 为 等 叶 高 直 叶 片 , 因 
此 在 展 向 上 流动 特征 相同 。 采 用 在 展 向 位 置 流动 特 
征 相同 . 图 6~8 分 别 为 展 向 为 50% 叶 高 处 的 流 线 图 
和 误 黏 性 、 压 力 云图 . 
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图 6 叶 高 50% 处 流 线 分 布 图 


Fig. 6 The streamline distribution at 50% span 
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图 7 叶 高 50%% 涡 务 性 云图 
Fig. 7 The Eddy Viscosity distribution at 50% span 


Pressure 


图 8” 叶 高 50% 压 力 云图 
Fig. 8 The pressure distribution at 50% span 


图 6 为 叶 栅 50% 叶 高 处 的 流 线 图 , 由 图 可 知 , 流 
道 通畅 , 在 静 叶 尾 缘 有 流动 尾 迹 , 但 没有 流动 分 离 产 
生 . 因 透 平 设计 压 比 在 0.14, 整 级 灼 降 较 大 , 动 叶 为 


缩放 形式 ， 在 动 叶 内 为 超 音速 流动 . 图 7 为 涡 昔 性 
云图 , 在 静 叶 出 口 调 务 性 有 少许 增 大 . 图 8 为 压力 云 
图 , 在 整个 流 道内 顺 压 流动 。 

3 变 工 况 性 能 预测 

3.1 进 汽 压力 变化 对 性 能 的 影响 ( 背 压 不 变 ) 

透 平 热源 的 变化 会 引起 蒸发 压力 变化 ， 当 进 汽 
压力 增 大 时 , 膨胀 比 增 大 , 级 的 降 增 大 . 在 设计 转 
速 和 背 压 下 , 研究 进 汽 压力 变化 对 级 的 性 能 的 影响 ， 
进 汽 压力 0.33~0.88 MPa,， 间隔 选取 0.11 MPa。 效 
率 、 流 量 及 功率 随 进 气压 力 变化 规律 如 图 9、10 所 
示 ， 进 汽 压 力 增 大 理想 灼 降 增 大 ， 效率 随 进 汽 压 力 
增 大 先 增 大 后 减 小 , 进 汽 压力 为 0.66 MPa， 即 灼 降 
30 kJ/kg 时 效率 最 高 . 在 进 汽 压力 变化 时 ， 喷 嘴 达 
到 临界 工 况 ,流量 随 着 进 汽 压 力 增 大 线性 增 大 ， 输 
出 功率 随 着 进 汽 压力 压 增 大 而 增 大 ， 


03 04 05 06 07 08 09 
Po/MPa 


图 9 效率 随 进 汽 压力 的 变化 曲线 


Fig. 9 Efficiency of vs. inlet pressure 


3.2 进 汽 温度 变化 对 性 能 的 影响 ( 进 汽 压力 、 青 压 

均 不 变 ) 

热源 的 温度 变化 会 引起 进 汽 温 度 发 生变 化 ， 当 
进 汽 温度 发 生变 化 时 , 级 的 烩 降 、 效 率 、 轴 功 也 会 相 
应 发 生变 化 。 进 汽 压 力 为 设计 压力 ， 当 进 汽 温度 从 
373 K 增加 到 408 K 时 , 间隔 选取 5 K, 工 质 由 饱和 
进 汽 逐 渐 到 过 热 的 过 程 ， 研究 进 汽 温度 变化 对 透 平 
性 能 的 影响 . 结果 如 图 11、12 所 示 . 

随 着 进 汽 温度 的 增加 , 效率 逐渐 下 降 , 但 是 变化 
小 于 1%。 随 着 进 汽 温度 增加 , 级 理想 灼 降 增 大 , 在 
进 汽 压 力 和 背 压 不 变 的 情况 ， 比 容 增 大 ， 质量 流 量 
减 小 , 但 由 于 灼 降 增 大 , 效率 基本 不 变 , 输出 总 功率 
仍然 增 大 . 因此 可 认为 ,离心 透 平 与 轴 流 和 向 心 透 
平 热 力 特性 相同 ， 进 汽 温度 变化 对 运行 经 济 性 影响 
较 小 。 
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图 10 功率 和 流量 随 进 汽 压 力 变 化 曲线 


Fig. 10 Output power and mass flow rate of inlet pressure 
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图 11 效率 随 进 汽 温度 的 变化 曲线 
Fig. 11 Efficiency of vs. inlet temperature 
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图 12 功率 与 流量 随 进 汽 温度 的 变化 曲线 
Fig. 12 Output power and mass flow rate of vs. inlet 


temperature 


3.3 背 压 变化 对 性 能 的 影响 ( 进 汽 温度 、 压 力 不 变 ) 

透 平 出 口 背 压 为 冷凝 温度 对 应 的 饱和 压力 。 当 
环境 水 温 随 季节 发 生变 化 、 循 环 水 泵 的 变 工 况 运行 
等 均 会 引起 冷凝 温度 发 生 改 变 ， 进 而 使 透 平 出 口 背 
压 偏 离 设 计 值 . 原 设计 背 压 0.11 MPa 对 应 的 冷凝 温 
度 为 30"C. 当 冷 凝 温度 从 27°C 到 45°C 变化 时 , 背 
压 从 0.101235 MPa 升 高 到 0.1822 MPa, 背 压 变化 对 
透 平 性 能 的 影响 如 图 13、14 所 示 , 效率 随 着 背 压 升 
高 先 增 大 后 减 小 , 在 背 压 为 0.134 MPa( 冷 凝 温度 为 


35°C) 时 , 即 灼 降 为 30 kJ/kg 时 , 效率 最 高 . 由 于 在 
该 背 压 变化 范围 内 ,， 均 在 临界 工 况 下 运行 ， 所 以 流 
量 保 持 不 变 , 但 随 着 冷凝 温度 的 升 高 , 背 压 升 高 , 理 
想 炊 降 减 小 , 输出 总 功率 减 小 ， 
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图 13 效率 随 背 压 的 变化 曲线 


Fig. 13 Efficiency of vs. exhaust back pressure 
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图 14 功率 与 流量 随 背 压 的 变化 曲线 
Fig. 14 Output power and mass flow rate of vs. exhaust 
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图 15 效率 随 转速 的 变化 曲线 


Fig. 15 Efficiency of vs. rotation speed 

3.4 转速 的 变化 对 性 能 的 影响 (设计 参数 ) 
在 透 平 运行 时 ， 会 发 生 转 速 偏离 设计 转速 的 工 
况 , 因此 透 平 要 求 在 一 定 转速 波动 范围 内 , 维持 良好 
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的 性 能 . 在 设计 热力 参数 下 ,对 转速 在 6000~15000 
r/min 范围 内 变化 进行 了 数值 模拟 , 效率 随 着 转速 变 
化 先 增 大 后 减 小 ,在 设计 转速 变动 25% 的 变动 范围 
内 ， 效 率 均 在 80% 以 上 , 规律 如 图 15、16 所 示 . 在 
热力 参数 保持 不 变 , 转速 改变 , 质量 流量 保持 不 变 ， 
输出 总 功率 随 着 效率 先 增 大 后 减 小 . 
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图 16 功率 与 流量 随 转速 的 变化 曲线 


Fig. 16 Output power and mass fow rate of vs. rotation 


speed 


4 结 论 


本 文 以 有 机 工 质 R123 为 流动 工 质 , 借鉴 轴 流 透 
平 气 动 设计 方法 , 设计 了 应 用 于 低温 热源 利用 的 离 
心 透 平 , 用 商业 通用 软件 CFX 对 所 设计 的 离心 透 平 
整 级 性 能 进行 数值 模拟 和 变 工 况 性 能 研究 . 

研究 表明 : 

1) 在 设计 工 况 下 , 流 道 通畅 , 顺 压 流动 , 无 流动 
分 离 , 效率 和 功率 达到 设计 要 求 . 

2) 在 变 工 况 条 件 下 , 当 进 汽 压力 增 大 , 灼 降 随 之 
增 大 时 , 效率 先 增 大 后 减 小 , 膨胀 比 在 5~6 之 间 , 也 
即 灼 降 在 30 kJ/kg 左右 , 效率 最 高 . 随 着 进 汽 压力 增 
大 , 流量 增 大 , 输出 功率 增 大 ; 冷凝 温度 在 27~45°C 
范围 变化 时 , 也 即 背 压 在 0.10~0.18 MPa 变化 时 , 级 
在 临界 工 况 下 工作 , 随 着 背 压 增 大 , 流 基 不 变 , 效率 
先 增 大 后 减 小 , 效率 仍 在 灼 降 为 30 kJ/kg 时 为 最 佳 ， 
输出 功率 随 着 背 压 的 增 大 而 减 小 。 

3) 进 汽 温度 改变 对 气动 性 能 的 影响 较 小 ,效率 
变化 范围 在 1% 之 内 , 功率 随 着 进 汽 温度 的 增 大 而 增 
加 ; 转速 改变 时 , 在 设计 点 效率 最 高 , 转速 改变 对 流 
量 没 有 影响 。 
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